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PRSFACE 

\ linear automatic control system must be stable and yet meet certain 
speci fications of accuracy. In veneral, voth transient and frequency 
resoonse are needed to guarantee satisfactory designs. 

The purpose of this investigation was to show that phase margin can 
be used to correlate transient and frequency response characteristics, 
primarily the resonant frequency peak and the transient peak overshoot. 

Using the electric analog computer, numerous servo systems were simu- 
lated, The resonant frequency peak ad transient peak oversnoot were plotted 
versus phase margin for various system gain settings. From the plots approx- 
imate equations relating the resonant frequency peak, the transient peak 
oversnoot ani the phase mrgin were developed. Soluti ons to the equations 
are rererally accurate to within a small percentage of the correct values. 

For desi™m purposes, a more useful correlation of resonant frequency 
peak and trinsi2nt peak overshoot exists by plotting the two versus a 
dimensionless parameter i » if curves for a sufficient number of systems 
are grouped in families, accurate resonant peak and transient peak values 
are determined easily with the knowledge of the system gain amd the system 
time constants. 

This report has show that phase margin is the paraneter that corre- 
lates the resonant frequency peak md the transient peak overshoot. The 
correlation may be partially ambiguous if the system possesses two or 
more vhase margins of the same value. However, the direct correlation of 
resonant frequency peak and peak overshoot is «till valid. 

System gain and gain mrgin also display a regular variation with 
phase martin, But the variation is dependert on other system character— 


istics as well as phase margin. 


ii 











Phase margin can be a powerful tool to the desigrer of servo systems : 
Its use is not limited to certain systems but is applicable to systems of 
higher order, It is easy to use am gives reliable results. 

The author wishes to express his appreciation to Dr. G.J. Thaler, U.S. 
Naval Postgriduate School, for his assistance, suggestions, cooperation 


and encouragemerit in the preparation of this paper. 
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TABL® OF SYMBOLS AND ABBREVIATIONS 


(Listed in the order of their use in the text) 


The output or controlled variable. 

The input simal. 

Time or transient response. 

Radians per second, 

Frequency response. 

Phase margin, equal to 180 degrees minus the 
transfer function phase angle at the frequency 
at which the gain is zero decibels. 

The |-. | resonant freouency peak. 

Transient peak oershoot. 

The system transfer furr tion. 

Phase angle of KG (iw). 

System damping coefficient. 

System natural resonant frequency. 

System gain. 

Tine constant. 

Gain mrgin, [ 1 KG( jw) a at frequency at 
which Ang G(jw) = -180°. 

Resistor, 


Capacitor, 
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“CORRELATION GO TRANSIENT 
AND FREQUENCY RESPONSE 


‘Y USE OF PHASE MARGIN 
CHAPTSR I 


INTRODUCTION 


™n the desiyn ind adjustment of linear servo systems the designer 
nust obtain a system that will be stable and will meet certain speci- 
fications of iccuracy. In gereral, both transient and frequency 
resvonse studies are needed to guarantee satisfactory designs. 

The performance of a control system is frequently judged on the 
basis of its response to a transient input consisting of a step function 
of position. Abrupt disturbances represent important conditions under 
wnich tre system response is of interest and in which knowledge of 
svstem oerformance on the basis of sinusoidal ats alone is not suffic- 
jent. Because of the ease with which the system response to a transient 
innut nay be determined experimentally «ith limited equipment. it is 
desirable to be able to check design results against actual system per~ 
formance. 

Although dynamic stability requirements are more easily based on 
transients followin’ step and ramp disturbances than on the behavior 
following sinusoidal disturbances of varying frequency, the design of 
conventional systems to fulfill the majority of performance specifi- 
cations is more easily accomplished by studying the response to a 
sinusoidal disturbance of varying frequency. An important reason for 


the use of the sinusoidal procedure is the relative ease with which 





the dynamic properties of an intricate servomechanism can be ‘represented 
by an experimentally determined frequency function, while representation 
of all its properties in terms of a differential equation may be nearly 
impossible, 

The effectiveness of the sinusoidal procedure involves correlation 
between freouency functions am time functi ons. Correlation to a degree 
has been based on noting: 

(1) salient features of the transient such as speed of response, 

“requency of oscillations, degree of stability (i.e., peak 


overshoot), and steady state error, and 
Oo 
Se 


the phase angle relating error to output at low frequencies. 





(2) resonance peaking of | at certain frequencies, and 

The transiznt response of a system may be determined analytically 
oy the Laplace transformation. But an exact transient analysis for 
conpl3x systems involves an undue amount of labor. Frequent use is made 
of dynamic analyzers and analog computors to limit the time required for 
Ootaining a solution to transient am frequency responses. 

An appRewe: to determining the transient response of a system is to 
determine some correlation between the nature af the closed loop response 
of systems to steady state am transient inputs. Thus,the transient 
response of a system having a known response to a steady state input can 
be determined, 

If the frequency response F( jw) is known, the transient behavior 
f(t) in the time domain can be found. Conversely, if the time response 


f(t) to a unit step input is known, then the response F(jw) in the fre- 


quency domain can be found. The responses are related by the integral 





equations: foe) : ' 
| w r 
f(t) = 3n/_' é dw Ce 
and fore) -jwt 
Fijwi=/ fltie dt (2) 
— 00 


These equations do not permit general explicit solutions. Their relation- 
ship may be explained in the following manner. An arbitrary input signal 
exoressed as a function of time can be represented by a Fourier spectrum, 
The rrequency response curves of the total system indicate the extent to 
which each frequency component in the signal is transmitted b:’ the system, 
i1.e., the change it undergoes in vhase and amplitude. The frequency— 
anplitude curve representing the input signal multiplied by the amplitude 
transfer curve of the system gives the Peurter spectrum of the output 
cipnal. The phise function representing the input signal added to the 
phase function of the system gives the phase furction of the output signal. 
However, the output signal expressed as a time function may also be ex- 
pressed by the Fourier synthesis of the component frequencies. Thus, it 
is peas@nable expect correlation between the frequency response and the 
transient resvonse curves fa@ the total system. 

Since correlation of a kind does exist between the transient response 
following a step input and th sinusoidal response as given by the fre- 
quency functions, the frequency response or loci studies are highly in- 
portant in the design of servomechanisms, 

The primary purpose of this investigation is to show that there exists 
for various types of servo systems a regular variation of phase margin 


with the transient peak overshoot P.O. and the resonant frequency peak Mp. 





The data was compiled by computative, graphical and analog computer 
studies. 

The primary result is a series of non-dimensional curves for er 
transfer functions, these being usable directly for design or gain adjuste- 
mnt purposes. By the use of these curves the designer can tell for an 
anpropriate system whether he can secure adequate stahility am still 
satisfy other specifications such as accuracy. Since steady-state accuracy 
and transient response are wholly interdepvendent and fixed once the loop 
gain function is established, the designer knows immediately if compensation 
is required. 

A secondary prupose of the investigation is to determine an optimum 
use of phase margin for higher order systems. Present practice generally 
involves use of M-Y contour charts [1] to predict resonant frequency 
peaks. The phase tmarein correlation should permit a relatively accurate 
prediction without resort to the M-N contour charts. 

A third objective is to jetermine a procedure for future experimental 
study to further develop the phase margin correlation such that all the 


a 


system requirenents of interest may be determined readily from a knowledge 
~*~ 


cd 


of the phase rargin. The correlation examined here involves only peak 
overshoot and resonant frequency peak. Other important characteristics 
include speed of response, frequency of transient oscillations and steady 


state error, 





CHAPTSR II 
PREVIOUS INVE&STIGATICNS 


Uxisting correlations of the transient and frequent resnonse of 

linear control systems may be classified as follows: 

(1) OCualitative, due to Bode [2], Ferrell [3] , Graham [4], 

Harris [5] and Hall [6] . 

(2) Analytic, for second and third order systems [7,8]. 

(3) Guantitative, due to Floyd [7] , Dawson [9] , and Harris [to]. 

(4) “mpirical, due to Chestnut and Mayer [1/J , Herr and Gerst [12] 

and Hunnicutt [13]. 

A brief discussion of the various methods follows: 

1. (a) The criteria of servo system stability utilized in the fre-— 
quency approached is based on H. Nycuist's [14] theoretical 
analysis of regeneration in feedback amplifiers. In the stabil- 
ity adjustment of a servo system it was found necessary to provide 
margins of safety because of aging and other factors. H. W. Bode 
(2] discussed the theoretical considerations involved in the 
selecti dh of such margins in the desig of feedback amplifiers 
and introduced the concept of phase and gain margins used extens— 
ively by communication engineers. 

Phase margin OM is frequently defined as the value of the angle 

[160° ~Ang G(jw)] for the frequency w at which KG(jw)=/ . Gain margin 

is defined as the log modulus of [/-Kk&G jw) wien the frequency w is such 
that the Ang G(jw)=-/80 . Neither aN nor gain margin is related directly 
to M,. Sinilarly, the mere existence of aQor, say, 40° does not insure 


p 


a satisfactory response. 





rT 


(bo) 3.3. Ferrell [3], in discussing the applic.tion of Bode's 
treatment to servo systems, stated that a phase margin of between 
LO degrees and 60 degrees and a gain margin of 10 to 20 decibels 
constituted good desig practice. Furthermore, Ferrell showed 
that the phase shift varied with the slope of the gain curve such 
that 15 degrses of phase shift accompanied a one decibel change of 
gain per octave. 

(c) R.~. Graham [4], stated that if the gain margin is sufficient, 
the phase margin is usually the dominant factor in caetermining the 
amplitude of the initial overshoot. The gain margin needed ina 
particular design depends on the expected variability of the loop 
transmission, Radar tracking loops usually lmve gain mrgins of 
the order of 15 decibels or more because of tle large number of 
factors which may cause loop gain to vary. 

(d) H.Harris, Jr.,{5]has stated that as the magnitude of the fre- 


KIT 


quency function changes, the phase angle approaches — “37 





radians for increasing frequencies. If the loop transfer function 
magnitude decreases as + or faster for a region near unity gain, 
the phase angle will be near minus 180° or more and the system will 
be highly resonant. A good loop transfer function will not decrease 
in magnitude faster than about a until the megnitude is below 
about one. Then the resonance peak will be no greater than 13. 

(e) At. Hall 6] gave an approximate index of the transient response 
based on the comparison of frequency response curves and the trans-— 
ient response for the ae curves. The presence of peaks in the 


amplitude response was associated with the complex roots of the 


characteristic equation. The height of the peak relative to the 





flat, low frequency part of the curve is an index of the real part 
of the root,tending to increase as the real part decreases. The 
angular frequency at which the peak occurs is an index of the imagin- 
ary part of the root, tending to increase as the magnitude of the 
imaginary part increases. Thus, the frequency of the peak is an 
index of the frequency of the oscillatory component of the transient 
response and the height an index of its damping. In numerical terns, 
if the height of the peak relative to its response at zero fre- 
quency is limited to 1.4, then the damping ratio should lie between 
O.5 and 0.8, and the angular frequency of the peak should equal 
the frequency of oscillation to within twenty percent. 

2. (a) G.S. Brown and D.P. Campbell [7], ami others, have shown that 
a correlation exists between the maximum or resonant amplitude of 
the frequency response Mp and the peak overshoot P.O. for various 
values of the damping ratio §& . As & decreases, the value of 
Mp increases and P.O. increases. The magnitudes of the peaks of 
the amplitude response are measures of the damping of natural modes 
of ori dation of the transient response. The frequency band over 
which the amplitude response has a substantially constant magnitude 
is a measure of the speed of response to transients since a high 

- natural frequency, and therefore a high speed of response, is 
linked with a high resonant frequency in the amplitude of the 
Sinusoidal response, 

For a system having a second order characteristic equation, an exact 


relationship for Mp as a function of § has been computed: 





- a : Ww (3) 
oo +j2s 





Differer.tiating ind setting ecual to zero: 


w= w, {1-257 (4) 








_| 6. ] _ | 
then Mp = S. |= SETH (5) 


It is 2asily shown that in terms of the same parameters: 
T 





i-ea 
PO, = |+ = sin (tan / 8 ) (6 ) 

\n examination of Mp shows that it becomes excessive for values of 
d04 isa result, widespread practice has been for ?.p to be limited to 
l2<Mp< 16, 

2. (vb) For the thiri orier system, no simvle relationship exists b2- 

teen system pirameters and the coefficients of the factored cubic. 
Charts relatine the coefficients of the general cubic 
dyx > +a,x*+a,x +1=0 (7) 
to the values of a,§ and w, in the factored cubic have been pre- 
pared bY @.J. Liu and L. J. “vans ° 
(Ho x + Ide + GPx +1) =0 ) 


a 
Tre use of the charts require appreciable. calculations since the 


coefficients ang a, and a, are not related in any simple way to actual 
System raraneters., Charts relating system parameters to system verform— 
ante nave been nade un by H.M. <smes ( 8]. However, no correlation of 
Po... ith the system oiram-ters vas mde, 
3. Analytic methods fcr correliting steady state and transient re- 
sponse bas2d on the Laplace or Fourier transform ‘eve bean de~ 
veloped. One method is that of G.*. Floyd described by 3rowm and 


Canpbell { 7). Floyd's method is based con an approximute solution 











of the inverse trinsform of the known frequency function. The 


nethod is easier to use tlan any exact mthod for systems greater 
than th: fourth order and is accurate to better than 5 per cent. 
More recently H. Harris, Jr., M.J. Kirby, ard 3.F. Van Arx (9], 
und C.H. Dawson [io] have published methods of developing transient 
response curves from frequency response. 


(2) Another method of determining transient resoonse makes use of 


= 


a series of chirts requiring no detailed calculations. These charts 
relate the open-loop response to either a sinusoidal input signal 
-r to 1 step input simal. 

H. Shestnut and RoW. Mayer [NJ have developed such charts for a multi- 
tule of systens of a cert2in character, mmely Type I with the third, 
fourti, or fifth order characteristic equation, 

4. (bo) D. Herr and I. Gerst[i2] have correlated the responses of a 

limited number of systems possessing minimum band-width. Herr und 
Gerst did not zraph their correlations but presented them in tubu- 
lar fein, 

ihe orien oe just described are limited in that they avply only 
to syste..(s ee orescribed open-loop attenuation characteristics and are 
not vereral. 

WR. Eunnicutt, Jr.fi3] measured gain margin, phase margin and trans- 
lent r2soonse for non-lin2zar systems subjected to amplifier saturation. 
Hunnicutt vlotted D,, versus various system parameters but did not plot 


= or P.9. versus G and did not arrive at any successful correlation. 











Ciarment LIT 


PROGEDUR 


Sy reference tO fame, 1, it is seen that wy is a fumtion of the 
corner fre uencies, tiie sloxe of the loop trinsfer frequency curve and 
tho ssten cain. fne loon transfer furtion uniquely uefines all of tie 
properti2s of the lin®ir servomechanism, novwever, for each suen transfer 
function, there is «n infinite variety of physical assemblies of servo 
estes havin this trimaier runction. Hence only the mathematical form 


of the trinsfer fumtion was considered here rather thin particular physical 


on.i Ur tions. 


Various s steis of th= second, third and fourth orier am of Tepes I 


joa TI werd lavectivtited with system Zain the narameter to be viried. Thi: 
is 2° Wivilert to th i~sien problem of adjusting systom zain to sive ie- 
mrei verforrance. For exampl?, tske the transfer fur tion: 
_ Ky (jw Tz +/) 
kK GG») (fe Get +t wT +) ©) 


moe COR MaC TSI AS scuntion is: 


. Kh kK (jo) 

aw tee b)Gu)*+ Get SE) jw +g = 0 
The syste 72in term thenges both the system darping ind natural frequency. 

& schematic dlasram of the equipment used to reisure the peak fre 
quency response ind the peak overshoot is shown in Fig. 2 am equipment 
Mmeelist-d in \ppemiix @. 

» loon tri.nsfer furtion wis first slotted on the Bode or decibel-log 
Wrlot, Various values of sycten rin “ere selected suwh tmt the phase 
rarvcin was varied from anout 2 iegrees to 70 ueprees. Mn wis determined 


by nlottine inoropri ite peints “n MN contour charts. P.O. owershoot wis 
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calculated for at least one point fa@ all systems examined that were of 
third order or less. 

The basic element in the system was the analog computer. A Boeing 
Analog Computer, BOAC, was available am this was used. Automatic balance 
units were not available but drift did not give trouble except for high 
gain systems. In these circumstances, voltage scaling was an adequate 
measure. 

The sine wave generatar used was the Mod.202A Hewlett-Packard Low Fre~ 
quency Generator. At low frequencies the irregularity of the wave shape 
was quite apparent but since the components contributing to this erratic 
wave shape were of a much higher frequency order, and very small in ampli- 
tude, they did not affect the response of tle low pass systems, 

The step input was supplied by 45 volt battery installed in the BOAC. 
The input and output signals were recorded by Brush Magnetic Pens. 

For the experimental tests, the R-C components used in the computer 
were IRC resistors, having + one per cent tolerance am a temperature 
coefficient of .025 to .065. Sprague Vitamin Q Capacitors and Cornell- 
Dublier Metal oe Dykanol paper capacitors were ued. 

The odnomtte were measured to four significant figures. It was as- 
sumed that their values while in use did not change. “xperimental values 
of P.O. for low gain systems not using voltage scaling differed from calcu- 
lated values by 0.1 to 2 percent. On higher gain systems deviations as 
high as 7 ver cent were recorded. 

Thre measurerent of Mp on the Brush recorder was not so good. Since 
the resonant mignification was frequently 3 times or mare, the initial gain 
settine of the Brush anplifier was so small tmt it could not be set and 
measured with sufficient accuracy. Hence, the values of Mp used in this 


report are thase determined on the M-N contour charts. 


Li 
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R4SULTS AND DISCUSSION 


The sivnificant data taken is recorded in TablesI and JI in Appendix 
3. Tris data has been vlotted in two different mamners: 

(1) Mp and P.O. versus @, , Figs,3, 5, 7, and 9, and 
(2) Mp and P.O. versus ty » a nondimensional parameter, 
Figg, 4, 6, 8,=and 10. 

The similerity of thevecurveenon the first type of plot may be quite 
3urprising. The sole mrked variation is in those systems in which phase 
lag or leid toxnpensation is inherent or has been intentionally added. 
This oc*urs whenever two different system gain settings can vroduce the 
sam-2 PQ, - Iven here, though, it is seen tliat a similarity exists in that 
the .isher sain setting always leads to a larger Mp and P.O. for the same 

e. .lso, ‘or the conditionally stable syste:s, it is aoparent that as 
Stavility and Q, decrease the curves of P.O. and Mp become indistinguish- 
able from those of systems that are absolutely stable. 


Ar aonroximate linear relationship between Z, and P.O. is: 


log, RO. = 0.259 (,- Ba) 0, 160 (11) 


7115 empirizal formul2 is accurate to better than lv percent for Drs near 
20 or -\ J2preaes and vives better than 5 percent accuracy for the region 
near 4) ie7vrees, Th looped back portion of the curves (increasing ge 
vith intreasin system zain) wis neglected in the preceeding statement but 


in most cises the empirical equation will still be accurate to within 10 


percert. 
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Similarly an aoproximite linear relationshiv between Q omer .O0, iS? 


log, Mp = = see) ,40°2 J 460° (12) 


ifain nefrleeting th folded back portion of the Z-Mp curves (inctreiwsing 
Z. with intrecsinz system gain), this emoirical equation is .iecurate to 
within <.5 percent at 40 uegrees and 5.0% at 6U dervrses vhase margin; 
tiese D,, Limi: s inclug@® the range of greatest interest. 

vith these two sinple equaticns a designer, knowing Q only, can pre- 
dict .ith batter than 1U percent aceuraicy the Mp and P.O. to be exp:cted. 
If tir Z, is tiken in the range where it is incr2asing with increasing 
system z1i1, these squations sive optimistic results, tie err-.r being 
greates as th 2 bind=-width of this characteristic incre3ses,. 

It is nossibdl2 to fit empirical equations more accurately to the mein 
Mp- Dn and >.°. - Oeieves out the improvement in accuracy doesn't warrant 
she aid tion.l complexity involv2d. Should a more exact exoression be 
desired, in inverse transcendental or a hyperbolical curve such as 

Z,, =A tanh (B log, Mp + C)+D (13) 


eS ould be tried. 
~ 


band 


The const wits A, 8, C and D are used to shape the curve [15]. 

In Fic. il, Mp versus P.O. is plotted. For tre systems exanined in 
this report, it is s2en that “p versus P.C. ioes have a general relation- 
st31D. main a lineur eupniricil formula is used to cezcribe this corres= 
ponae rice . 


4) 
logig Mp = 2.015 log, PO, — 0.087 





“or the swste s of this report, this equation is iccurat» to about 15 vercent. 


ewever, fori most systsms covered 1t is Wccurate to Within 1G percent. © 
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now have a simple relationshiv between Mp and P.O. such that if Mo or P.O. 
is ieteam.ined experimently, the other is known to within 10 or 15 percent. 

Having observed these correlations of ip and P.O. with @, and Mp with 
>.O., an analytic correlation would seem desirable and possible. Taking 
the expressions for ‘‘p and P.0O., equations (5) and (6) in Cnapter II, a 
relationship between P.O., Mp and 6 can be derived. However, the expression 
in the forn of a series in $ and ‘p does na& converge rapidly and so is of 
micbie practical aid, 

The question arises is to whether the empirical relations derived from 
Fire. 3 to 10 are useful for other and perhans more complex systems likely 
to be encountered. In Table III, Appendix B are listed data for the systems 
analyzed by Herr and Gerst [i2] . For example, taking a system listed under 
Case III, having a fifth orier differential characteristic equation: 

Z,=¢50° 5 Mpaia 5 RORLSS 

logy Mp = 0.435 (1- FEE) =0.15S 9, Mp 2156S 

% Error =-pe2 4/00 =2.50% 


login FO. 20.259 (1- Fe tt) = 0.// PP ora .2\9 





06 
% &rror = 735 x/00 =4h 4 5% 


If the vo is know, using the empirical relation for P.0O.: 
logig Mp = 2.018 X 13-087 = O17S »Mp=1F97 


.O9 
% Error = oe! x J00 =6.93 v4, 





To take an example of a fourth order equation listed under Case IV: 
OQ, = 60 Mp = Los Pe a a = | Fa 
GO ard 
loging Mp = 0.435 (1- Sez) = 0.0174 » 7p / 4- 
fay eds 
% Evror = Los X/00 = O. G5 % 
log), FO, ce 259 (i- <2 = 0.52. , PO.=1/28 


% Evror =-TF* 100 = 2.54% 
login Mp = 2.015 x,04] -.087 =-0,0045 9 Mp9 


06 
DECOY = aes x /00 =5,8/ We 
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By the use of two simple empirical formulas, values of Mp, P.O. 
or Z, cin be determined by knowing any one of the three. Usually these 
vilu2s ar? accurite to within 10 percent or better. 

Generally the designer is given a set of specifications from which 
he can determine the required steady state frequency response and the 
transient response to a step input. Frecuently certain elements of the 
servomechanism are fixed as part of the system and cannot be varied. Also, 
tre loid or driven element is often a fixed component of the system. With 
the fixed components plus a knowledge of servo system -elements, the de- 
Signer can determine the remimer of the componets necessary to obtain 
a suitable overn—-lLoop characteristic. It is a simple mtter to draw the 
ohase angle curve especially if a tangent template is available. Oms 
the time constants of the system are fixed, Mis changed only by chang- 
ing system eain. By the use of the empirical equations (11) and (12) the 
desicnrer has, with very little labor, a reasonable check as to whether 
system gain and dynamic stability are compatible for the system he has 
selected. If families of curves, such as Figs, 3, 5, 7, amd 9, are avail- 
able he has an accurate check on Mp and P.O. Furthermore, if families of 
mevyes, such as“ are indicated in Figs.4, 6, 8, and 10 are available, it is 
now Beet sayy to draw the phase angle curve, since in Figs.4, 6, 8, and 
10 Mp and P.O. ars plotted versus the dimensionless parameter i . Of 
course, the attenuation characteristic »of the transfer function must be 
Similar to one of the families since LY versus Mp or P.O. varies widely 
witn the virious op3n-loop characteristic curves. 

With the rinimum effort tne designer is able to tell avproximately 
by the use of empirical equations, or exactly by the use of families of 


curves, if specifications of stability and zccuracy are satisfied or if 
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counensition is required. The use of the empirical equations is general 


and not restricted to any particular system. The equations may be opti- 
mistic in the region previously noted, but that region is readily appar- 
ent from the open-loop attenuation characteristic. Tm error in the region 
of increasing phase mrgin with increasing gain is proportional to the 
band-wiith of the cornpensition. However, if either Mp or P.O. are know, 
tue other may ce conputed by equition (14) with the same degree of accuracy 
as stated before. 

Thus a regular viriation of Mo and P.O. «ith Q, does exist and 


may be used to correlate the steady-state and transient characteristics of 


4 systen, 
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“BAPTER V 


SUGG4STIONS FOR FURTHER 
INV«STIGATIONS 

In tr is riport ohise “airgin 4s bzen usei to correlat2 only two of 
m2 steagpestit] 44d ands proo2rtizs of a systean, Other prop=rtiss of 
interest are tine of resnonse to a step or ram») inout, time for trinsient 
to sett’: within prescribed error limits, the fresuency of transient oscil- 
latio., th: frecuency at which the resonant peak Mn oceurs, and steady 
meatle acewéhey, ~- 

Trat a remilar chinge in system gain accoxpanies 2 change in Or 
Js Shon in Fis, 12. However, the amplitude of the gain and its slope 
iivfars a creat daal with the virious systems. This is, in mrt, depend-nt 
on the voresence of a zero in the right half of the complex plane. The 
ritios of the tine constants or poles of the system also influence th? 
3ySt3an °lin versus Z,, relati onshin, 

In fice lg the ratio af gain margin divided by the lowest corner fre~- 
gueney <, is,plotted versus D, o Again a regular variation of gain 
Parcin «ith ones rarein is apparent; but like system zain, the gain margin 
sloca and “ai ~znitude devemis on other parameters as well as Dm ° 

The effect of the spacing of corner frequencies (or time constants) 
on D,, stould be further investigated. If the ratio is much greater than 
ten (or om cycle on the log w plot), the larger corner frequencies have 
little offect on g,, . Thie is equivalent to the correlation of Y. Chu 
[16] and J.G. “ruxal [17] showing that verfornance may be determined primar- 
il, by a complex vair of voles, simple poles near the orgin, uni the zeros 


in tne complex pline. 
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in this revort no systems incorporating varallel compensation in the 
reedbatk loop wer consiijered. Systeis having this type of compensation 
should be investivited to determine if the correlation of stability with 
Dy can be mide2 fully general. 

1S discussed in inpendix A, the experimental methad could be much 
inprovei, « dynwmit unalyzer is desired having constant amplitude step 
ani Ffreouency outputs with a frecuency range of O.1 to 1000 cycles per 
seconl, if t's analyzer meisures iccurately the phase shift and amplitude 
of the sinulited system the 2x» »2rimental technique is simplified and the 
accir cy of neisuremerts improved. In essence, then, a dynamic analyzer 
desiued to ooer.t2 over the optimum freouency range as determined by the 


ara loz connuter is needed, 








S— 


CHAPTER VI 


SONCLUS IONS 


In the desixyn and adjustment of servo systems, the goal of the de- 
simer is to develop a system that will be stable am will meet certain 
specificitions of accuracy. If the open-loop attenmtion characteristic 
is selected, the phase shift characteristic of the system is defined. To 
determine if phase mrgin, the suoplement of phase angle at zero decibels 
attenuation, could correlate dynamic and steady-state response was the 
purpose of this investivation. 

It has been shown that a regular am ordered variation of Mp and P.O. 
“i th Dy does exist. The correlation is displayed in Figs, 3 to .10 and 
an aoproximate linear relationship was show. The relationship was mde 
easier to use by plotting Mp and P.O. versus ~: » a dimensionless para- 
meter, suffering a loss in generality thereby. 

To the extent covered by this report, the correlation exists for 
‘systens having characteristic equations up to the fifth order amd should 
be equally COMmmSH. oT. for systems of all orders. Tle one ambiguity exists 
wena siven Z,, may exist with two or more different system gain settings. 
In this c1se ip is still related to P.0. by the equation derived by elimin- 
ation of J, from eountione (11) and (12). 

Future investigation must incorporate other system parameters to 
effectively correlate other dynamic and steady-state characteristics with 
phase margin. fA dynamic analyzer with a frequency range of 0.1 to 1V000 
cycles per second used with the 2lectric analog computer can simplify 


measurement techniaques,. 
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Corr=1.tion of st2zbility with phase margin should be of great assist- 


ance to the iesizmr in developing linear automatic control systems, 
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APPUNDIX A 


DXPARIMINTAL EQUIPMENT 


The experimental system is shown schematically in Fig. 2. The Boeing 
\naloz Computer BOAC was used to simlete the servo system. Automtic 
Balance Units were not available but drift gave difficulty only with higher 
gain syste7s in which voltage scaling was necessary. It is desirable that 
two persons work together when using the computer. Set-—ups are more quick- 
ly checked far mistakas, comouter checks performed am critical coeffic- 
ients determined. Theres is no substitute fa@ good computer technique. 
Peing pressed fa time, ao 2ffort was made to cut some comers. Partly 
as a result, an averace of three test runs were made for each good test 
run. ‘The troutles in mot faulty runs were detected, but in a few the 
cause of the difficulty was not found and the Semeteer set-up was rede- 
Sisred. in example wis inherent instability in a second order system et 
all zain settin,;s——but a second order system should be absolutely stable. 
Tre limits recomaended for system time constant components on the computer 
are 50K R& 20M and 0.000lWf£CZjouf. In gereral the range shaild be 
somewhat mor? Penis ota, Say, O/MERSIOM ard Moo meaCSuit . 

It was not possible to use the Mod. 1lO0U-A Dynamic Analyzer built by 
the Industrial Control Company since accurate readings of phase angle am 
magnitude tannot be mde at low servo system frequencies, 

The sine wives generator used was the Mod. 22-A Hewlett-Packard Low 
“remency Generator. The sine wave is formed by clipping a triangular wave. 
\t very lov frequencies the irregularities of tre wave shave were quite 
apparent, but the freouency camponmts of tr irregularities were of much 
higher order am s7zall in amplitude ae did not affect the response of low 


pass systems. 


BH 








The step input was supplied by a 45 volt battery installed in-the BOAC. 
If the battery voltage is low due to aging, it will be apparent on the brush 
recorder by abrupt random drifting in steady state. 

The input and output signals were amplified by Brush amplifiers and 
recorded on paper with Brush magnetic pens. Much of the work done was in 
the range of 0.1 to 10 cps. The frequency response of the Brush maghetic 
pens is about: 

0.2 cps—~8mm peak to peak, 

O.5 cps—lOmm peak to peak. 

80.0 cps—lOmm peak to peak. 

| In meisuring Mp, it is desirable that the recording instrument have a 

constint gain-frequency characteristic in the region of the resonant peak. 
The two pens used were not calibrated, tut it was assumed their character- 
istics were essentially the same. This was confirmed by the fact that no 
measurable differance in result occured irregardless of which pen was re- 
cording output ard the other input. 

The resistors used were IRC, with one percent tolorence and temper 
ature coefficient of 0.025 to 0.065. Sprague Vitamin @ Capacitors and 
Cornell-—Dublie Pigctal Cased Dykanol paper capacitors were used. The Vita- 
min @ capacitors were determined to hold a charge better than the Dykanol 
Capacitors. Polystyrene capacitors are much better than either of the 
above but were not available. The components were measured to four sig~ 
nificant figures, It was assumed that these valves did not change while 
in use, 

The reader interested in servo system instrumentation and masure- 

r mts is referred to F.C. Fickerisen and T.M. Stout [19], wW.S. Pritchett 
and R.if. Saunders [20], Rod. “hret (21J] or F.H. Ferguson and C.H. Looney 


[22] . F.B. Anderson [23] described a ten cycle to 10 MC analyzer. To 
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work «with the analogue c Ler a dynamic analyzer for measuring gain 
and nhase anvle through the frequency range of 0.1 to 1000 cycles is de- 


<a ravle, 
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APPSNDIX B 
“VELOPMENT OF DATA 

The systems tested in this repa@rt are listed and classified in Table 
Tj In Table 11 neem @ucr 2, Mp, P.o., 8 and ate 
as determined by tests, plotting and computation. The mthod of plotting 
was stated in Chapter III am is common practice. Tle calculations, where 
used, were also routine, albeit tedius. The plotting of Mp and P.O. versus 
@,, is apparently unique. Mp and P.O. vary in an exponential manner with 
D,, . However, in the region of greatest interest (40°S CaGO” ) ’ 
the curves may be avproximated by vem ient lines, equations (11) and (12), 
introducing only small error. The curves are shown in Chapter IV, pages 16 
Ee 23. : 

\{ more convenisnt method of using the data measured in this investi- 


gation is to plot Ip and P.O. versus “ o Knowing only the corner fre- 





quencies of the amplitude cmracteristic, system gain may be selected to 
satisfy dynamic and accuricy requirements. If this is not possible, the 
designer knows at the very start of the design problem that further compen- 


x 


sation is necessary. These curves are show in Figg, 3 to 10 in Chapter IV, 
my 


os 


pages 16 to 23. * 

In Table III are listed 4» Mp, P.O., and G.M. for the higher order 
sys tems. Mp and P.O. versus oO; are plotted in Fig. 14. Fig. 14 differs 
from Figs, 3 to 10 in that not only is ory varied but also the rutio and 
number of time constants. 

In Fie. ll Mp is plotted versus P.O. It is apparent that Mp has an 
ordered variation with P.O. The correlation is through their relation- 
ship with Z., » "saQuation (14) was derived merely by eliminating Dan 


in equations (11) and (12). This procedure turned out to be more accurate 
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than any relationshin that could be developed from Fig. 1l. 

The ratio of system gain to the lowest carner frequency of several 
representative systems is plotted versus D,, in Fig. 12. It is apparent 
that ths manner in which th system gain varies with Dna cannot be so 
easily correlated as were Mp am P.O. 

The ratio of gain margin to the lowest corner freqiwency is plotted 
versus Dr for several conditionally stable systems in Fig. 13. Like in 
Fig. 12, a correspondence is apparent but it must vary significantly with 


other parameters as well as m ° 
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TABLE JI 
SaaS we JESTEL 
Second Order Character/stic £9 ua #ion 





Ident. No. Ty pe Zero 
_ Kv aan 
mo ROG) = Cat eKjw ht!) wa! 
2-2 Same $70 
Ky ([ 4 7% +/ oe 2 
‘ald KGG4) *ConeDgwe ntl) sm! Bees 
2-4 Same ES, He/0 
ae S) Sam e A — $50 
Type One 
coe 
26 KEG) = GoGurel) , 
__ kv (jwia tl) di 
an? KG (jw) JHC wa +l) A aa 
2-8 Same Fr z/0 


1,3 





te 


(Aa a (Continued) 


Third Oreee Characteristic Lua tions 


74 Pe Zero 
Ke Gt ‘w/) “Gu ae (jw itl X { /w +h) 


Same 


Same 


/y pe One 
mae kv 
KG Go) = Gus Gn T4001) 


Same 
Same 


Same 


Ae bith 
KG Gu) = Je ett iju +/) 


Same 


Fourth Order Characterishic LE: 
Ty p@ Ore 


een ky Cw itee/) _ 
MG? GG FV GeRr 9% 


es kv Gw te +l) 


Lh 


guations 


Ss ,#+/00 


ay 2 Ts 
77, =83 *7, ~100 





System 


2-/ 


2-35 


Ge 
SAS8 
.28/ 
398 
FO 
.63 
sO 
/. 58 


57 O 
1.35 
10.0 
15:8 
an.2 
3/6 
50. 0 


0.25 


/.0 
58 
25° 
OX, 

1/00 
IEE 

25,0 

500 

100, 0 


PAE fi 


45 


Ves) LATA 
Dm deg Mp PF. O. 
67.5 = L417 
49.5 //9 128 
4.0.57 | 4F- 1.33 
36.0 1,58 137 
32.0 (8 /4/ 
26.0 2-24 LATS 
20.0 2.98 | $7 
59.0 1.02. 1/44 
ATS fae 28 
92,5 /.38 A276 
335° 176 / 366 
26.5° 2,2.4- L475 
24.0 J3.5/ L497 
19.0 S./2: 4576 
$65 07 40003 
43.0 /.38 VOT 
350 168 /238 
3S 488 L336 
30. O Zo [3597 
3/5 /. O4- 14/2. 
37.0 178 139 
435.0 ASF 3S” 
$1. S 1.46 AWB 
oF: 0 126 A278 
75.0 | lo 4224- 


Ass. 
Stab/e 


Abs. 
Stable 


Abs. 


Stable 


System 


2-4 


2-3 


2-7 


3.98 & 
/0.0 yo 
20.0 75.0 
$2,/ 75.0 

/00,0 17-5" 
3/6.0 82.5- 

12.0 52,0 
VG 39.0 

50. | 35S 
109.0 36.0 

25/0 44,57 
50/,0 S50 

S65 63.0 
/.0 57.0 
L76 40.0 
3./6 30,0 
I6/ 22.5 

/0.0 174.5 

10 540 

158 46.0 

2.57 39.0 

3.98 340 

6.3 30.0 
10.0 28.5 
15-8 28,0 
25,/ 29.0 
358 30, 0 


Mp 

0.978 
S94 
L005 
| O34- 
| O3F- 
/022. 


1.147 
153 
168 
1.72 
153 
4,36 


1,0 


LES 


h 436 
G6 
2.66 
3.6/ 


LMS 
[29 
1 SZ 
| 74- 
AIS 
2/1 
2./8 
2./4- 
2.02. 


1,6 


Fre £ LL (continued) 
Ky, Dm ,deg: 
7 


P.O. 
(,9004 
/.034- 
1055 
/,0S2. 


GM, 
WW, 
Abs. 

Stable 


ite 


Abs. 
Stable 


1/66 
278 
1.32.4- 
4328 
(1286 
[228 


;OSS Abs. 
1162. Stable 
1279 

1398 

|.507 

1 6os-— 


4237 Abs. 
129 @©Stable 
55 

A4/ 

/. #6 

1S 

/ $3 

|, $42. 
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System 
2-7 (cont) 


rire | 


3-/ 


3-3 


_s 


ti, Dy, de 


63.0 
/00.0 


/58.0 
250.5 


1-0 
2.0 
5,0/ 
/0.0 
25./6 
58 / 


1/90. 0 
200.0 
3980 


TSE 


/2.6 
20,0 
5/6 
SO, / 


6.3/ 
/0.0 
15.8 
25. / 


2.51 
3.98 
50/ 
6.3/ 
/0.0 


34:0 
38.5 
44.0 
51.0 


S70 
45-0 
36.0 
340 
I7-0 
45:0 
SIS 
67S5- 
77.0 


60.0 
47.0 
55, O 
23.0 
/3.0 


52.0 
38.0 
25:0 


15-0 


70.0 
420 
35:0 
3/.0 
19.0 


Mp 


188 
| 74 
1357 
/45 


4/0 
|. 34- 
/,64- 
18 
74 
[353 
1,35 
$22 
113F- 


L.Of 
[26 
166 
2.5/ 


4.47 


LIS 
IES 
2.37 
3.8 


ag 
1.26 
1.66 


1885 


3.09 


47 


F:@. 
1518 
1483 
1.441 
L383 


1129 
4232 
A337 
1.37 

1,356 
f 3 2. 
A225 
1177 


1/45 
4256 
1386 
S36 
4677 


/:2.4- 
36 
147 
A469 


[2.2 
4 356 
A438 
|. 504- 
£63 


TABLE ID Ceontinued) 


A bs. 
Stable 


A bs. 
Stable 


3.94 
reed 
1,73 
L#l 
AAS 


2-64. 
/. 84- 
L4-7 
(26 


2.6/ 
18/ 
62 
L465 
4272 


Sys fem 
3-4 


WA, Dn, deg: 


0. 56/ 
0.635 
0.793 
Vise) 
+26 
77 
2.0 
3./6 
4:/4 


0.3598 
0.63 
40 
A158 
Aon 


0.3/6 
OS” 
0.707 
1.0 
458 


0.3/6 
0.598 
OSs 
O56/ 
0.63 
0,75” 
10 


6/.0 
58.0 
52.0 
47,0 
4-/,.0 
34.57 
26.0 
2140 
15.9 


66.0 
SS" O 
42.5 
30,0 
19-0 


66.0 
550 
4505 
34.55 
205° 


So 
52.0 
45-0 
40.0 
36.0 
30,0 
20.0 


Mp 
4022 
Lost 
LIS” 
[242 
1.4 / 


1716 


2.09 
2.65~ 
3.76 


10 
/ 
138 
2.0 
3.02 


[0 
h/ 
AS 
470 

3.02 


LOS 
Vist 
{33 
[51 
47 
1,99 
2.98 


Po. 
LO6 
{1/0 
hls” 
Jef 
12.6 
1-36 
140 
1, S5~ 
463 


1,002.2 


O93 
123 


L39 
156 


103 
{MM 

jes 
136 
| S4- 


110 
L163 
125 
[29 
133 
143 
SST 


TABLE ID Ceontinaed) 


G. i, 
ASO 
7.56 
6.00 
SOC] 
3.80 
2.68 
| 74- 
156 
L2gs 


52.8 
2.86 
19/6 
4$2 
[1:30 


3./6 
2.07 
67 
|. 4-4- 
{26 


2.067 
A768 
4 $8 
[SOS 
|. 44- 
A315 
1426 








System ae, Pt, de ? 

3-8 40 55.0 

2,0 44.0 

5. 0/ 35S 

/2.6 F4.0 

257 $3570 

50./ 40.0 

100.0 49.0 

[$58.0 $53.0 

257.2. 56.0 

3I8BO 56.0 

63/.0 IO, § 

{0000 44,57 
2000.0 335.57 
3/60.0 27.0 

3-9 = /58 53.0 
3/6 $4.0 

Po 9G5/  2—_40.0 

/2.6 40.0 

25°/ 42.0 

SOS 4/.0 

/00.0 3570 
200.0 Zio 
4-/ 10.0 4-8,5° 
/00:0 SIO 

3/6.0 IF,0 

¢--2 159.0 0.0 
/000,0 32.9 
79000,0 37.5 


TABLE I (con tinued) 


Mp Ra SF, 
AJL 6 6ffas Abs. 
1355 {~ 254 Stable 
(66 £352. 
/80 /37 
(80 37 
/62, /3#F6 
44 {288 
43S ~252 
£27 (2lF 
Jam #17 
/20 
AIS 
172. 
2.// 
1/98 sl&# Abs. 
1,38 7/257 Stable 
[St 296 
15/5 432 
Sf L272 
LISS 224- 
466 £3523 
24/4 f#29 
[2 118 Fas 
[25 f23 TGS 
/77 f° 36 2, $2. 
Unstable 
49 L 45° 
‘Lee SE 


49 








£E/- OF Ff EGE ee gey sy &bE 


[p/-~ S584 f/ OCP OFZ pey 7 1h 
Ep/- O27 sY Gr Sisz. ofyY FY £9 
G9/-SE2/ 2/ tr E92 gz7/ 2/7 Fos 
Gia s/¥ /7% 929 wse Lt/y // #28 
AL/- C1) S07 Ie Lge bof s0/ IbS 


6z/- 907 0/7 9'29 GLS pyre! Jovy &z9 


a 4 ‘y /) ~~ 3 

x LP ade Ole 2(™!) = Pay alt LM; a> Za , 

OT aay aM I8N Oe eg? = (my Il 95°9 
3 


LE sy tbe 2 LLE/V $/ 4bE 


LEY pl Sep ' 292 #857 F/ st 
G7/ CY Yop 9°92 LEZ/ SY 99¢ 
fa) Sty. Jf £62 Lezy 2/7 wtys vA 
SHY {4 tL LES 8// fy OLS 

74s L60/ s0/ O19 £8f E60/ Soy ty 

sgy tof /0/ _S$9 2zas £50/ /o/ 949 

Fun on Ww Mpg PWD OL WW AP 


(uml mh fh < 
Fmt) Ay =( DOY . j®# a) 


C27 SWILSAS HLOIM-GNEOG WhaAW/NIWV 
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